Phase Diagram of the Two-Leg Kondo Ladder 
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The phase diagram of the two-leg Kondo ladder is investigated using computational techniques. 
Ferromagnetism is present, but only at small conduction electron densities and robust Kondo cou- 
pling J. For densities n > 0.4 and any Kondo coupling, a paramagnetic phase is found. We also 
observed spin dimerization at densities n=l/4 and n=l/2. The spin structure factor at small J 
peaks at q=(2n, Q)ir for n < 0.5, and at q=(n, l)ir for n > 0.5. The charge structure factor suggests 
that electrons behave as free particles with spin-1/2 (spin-0) for small (large) J. 

PACS numbers: 75.30.Mb, 71.10.Pm, 75.10.-b 



Numerical studies provide valuable unbiased informa- 
tion about strongly correlated electronic systems. How- 
ever, current computer limitations restrict most inves- 
tigations to one-dimensional (ID) models or small two- 
dimensional (2D) clusters. Unfortunately, the physics 
found in ID models is often qualitatively different from 
results observed in real materials. A possible procedure 
to start an investigation of realistic 2D models is by cou- 
pling ID systems together. This approach has been used 
with great success in the case of the t-J model, where 
numerical methods such as Exact Diagonalizationi and 
Density Matrix Renormalization Group (DMRG) 2 - were 
helpful in elucidating several properties of the 2D model 
through studies of t-J laddersi 1 ^' 4 ^ 

The ground-state properties of the 2D Kondo lattice 
model (KLM) remain mainly unexplored using unbiased 
methods. The KLM is the simplest model believed to de- 
scribe heavy-fermion materials^ and, thus, a better un- 
derstanding of its ground state is much needed. The 
main goal of the present work is to provide the first steps 
toward determining the phase diagram of the two-leg 
Kondo ladder (2-LKL). The information provided here 
will hopefully be as relevant for KLM 2D systems as lad- 
der studies were in the t-J model context. The TV-leg 
Kondo ladders consist of N Kondo chains coupled by the 
hopping term. As shown below, the phase diagram of the 
2-LKL is fairly different from the ID KLM. We consid- 
ered the 2-LKL on 2x1 clusters and a Hamiltonian 



H = -t E ( c l )CT s>+ H - c -)+ J £ s ^ 



3 ./ 1 



<i,j>,n 



where c J(T annihilates a conduction electron at site j 
with spin projection a, Sj is a localized spin-i oper- 



xj3 C j,a (T af3 



c- ~ is the conduction electron 



ator, Sj=|E ai 

spin density operator, and cr a p are Pauli matrices. Here 
(ij) denotes nearest-neighbor sites, J>0 is the Kondo 
coupling constant between conduction electrons and lo- 
cal moments, and t=l fixes the energy scale. The to- 
tal number of conduction electrons is N and n=N/2L. 
This model was investigated with the DMRG technique^ 
using open boundary conditions. The finite-size algo- 
rithm for sizes up to 2xL=80 was applied, keeping up 



to m=1200 states per block. The discarded weight was 
typically about 10~ 5 -1CP 9 in the final sweep. 

Let us first briefly describe what is currently known 
about the KLM ground-state phase diagram. In ID, 
for low electronic density and/or large J the ground 
state is ferromagnetic!^ (see also Ref . 0) ■ The rest of 
the phase diagram is characterized by a paramagnetic 
phase, except for a small wedge of ferromagnetism for 
fillings above n=0.5i^ Furthermore, spin dimerization has 
recently been discovered at n=0.5p4 Most ground-state 
investigations in higher dimensions have been limited to 
approximate approaches. DoniachJ"A> pointed out the pos- 
sible existence of a Kondo-lattice quantum critical point 
(QCP) due to the competition between the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction, which favors 
antiferromagnetism (AFM), and the Kondo effect, which 
favors paramagnetism. The full mean-field phase dia- 
gram of the 3D Kondo lattice was obtained by Lacroix 
and Cyrotpi 2 - They found that, at small J, there is a crit- 
ical density n c separating a ferromagnetic phase from an 
antiferromagnetic one. For sufficiently large J, however, 
the Kondo effect dominates and the system is param- 
agnetic. Further studies also considered an explicit ex- 
change interaction between localized spins^ Recently, 
quantum Monte Carlo (QMC)^ and DMRG±£ investi- 
gations of the half-filled Kondo lattice in small clusters 
confirmed the existence of a QCP at J~1.45 in agree- 
ment with previous approximate approaches*^ More- 
over, DMRG results on TV-leg Kondo ladders at half- 
filling have shown that the spin and charge gaps are 
nonzero for any number of legs and coupling JM> A two- 
channel version has also been studied at half-fillingiA£ 
Here, the two-leg Kondo ladder away from half-filling is 
considered. 

In Fig.^a), the schematic phase diagram of the 2-LKL 
is presented. We have identified three phases charac- 
terized by full ferromagnetism (FM), partially saturated 
FM (PFM), and paramagnetism (PM). The approximate 
boundaries between these phases were first obtained from 
the energy difference AE=E((2L - N)/2) - E(0) for 
L = 16, where E(p) is the ground state energy in the sec- 
tor with total spin projection <S'f = Si + sf=p. This 
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Figure 2: The spin structure factor S (q) vs. q x for the 2- 
LKL: (a) J=0.8, n=l/4 and L=32; (b) J=0.8, n=7/8 and 
L=32; (c) S (q x ,q y = 0) vs. q x for several values of L with 
J=0.8 and n=l/4; (d) S (q) for densities n=7/16 and n=7/8 
(see arrows) with L=32 and J=60. The solid (dotted) lines 
correspond to q y =0 (q y =ir). 



Figure 1: (a) Phase diagram of the 2-LKL. FM, PM, and 
PFM denote regions with ferromagnetism, paramagnetism, 
and partial ferromagnetism, respectively. To the right of the 
dotted line, the thick lines separate three regions: small and 
large J and n ^ 0.5. These regions are characterized by 
the location of the spin-spin structure factor peak (see text), 
(b) Gap AE between the ground state and the ferromagnetic 
state vs. J for a 2x16 cluster, and several numbers of elec- 
trons (as shown). The densities from the top are, n=0.1875, 
0.25, 0.3125, 0.374, and 0.5625. (c) Magnetization density vs. 
Kondo coupling, at several densities n. The errors are of the 
order of or smaller than the symbol size. 

was then double-checked through the (computationally 
more costly) calculation of the ground-state total spin of 
2x8, 2x16, and 2x32 clusters. In Fig. QJb), we show 
AE as a function of J for the 2x16 cluster, and several 
values of N. These values suggest that for small den- 
sity and large (small) J the ground state is FM (PM)jl& 
while for the whole region with density n > 0.4 and 
J>0 it is PM. In Fig. QJc), the magnetization density 
to=S't/2L vs. J is shown for some conduction electron 
densities. It can be seen that for small n and J the total 
spin is zero (within the DMRG precision). At n=l/8, 
m starts to increase as J is increased, and saturates at 
m=(l — n) /2. For n=l/4, however, m is non-monotonic 
and vanishes with further increase of J, apparently con- 
tinuously. These results are not due to finite-size effects: 
the same behavior is observed for both 2x16 and 2x32 
clusters at this density, suggesting that it survives the 
thermodynamic limit. For 0.25 < n < 0.4, the magne- 
tization density does not saturate at m=(l — n) /2 (see 



n=0.3125, Fig.Q](c)). Then, this phase has partial ferro- 
magnetic (PFM) order. For n > 0.4, we have found that 
m<0.03 for several Kondo couplings J (while for large J, 
m<10~ 3 ). This result strongly suggests that the whole 
region with n > 0.4 is paramagnetic. This is different 
from the ID KLM, which shows full FM at any density 
for large enough It is surprising that coupling just 
one more chain to the ID Kondo lattice induces such a 
dramatic effect on the phase diagram. 

To probe the paramagnetic phase, we calcu- 
lated the Fourier transform of the spin-spin corre- 
lation function (the spin structure factor) S (q) = 
2T E*,* e***-* <S£ • S£>, where S£=S ft +s ft . We 
observed that the maximum of S (q) can appear in three 
distinct positions, as indicated in Fig. [Ha). F° r small val- 
ues of J, the maximum of S (q) is located at <f=(2n, 0)7r 
for n < 0.5, while for n > 0.5 it is at q={n, As 
an example, in Figs.EJa) and (b) the spin structure fac- 
tor S (q) is presented for the 2x32 cluster with J=0.8 
and densities n=l/4 and n=7/8, respectively. These re- 
sults do not seem to be caused by finite-size effects, as 
shown in Fig. [2c): the peak becomes more pronounced 
as the length L increases. Previous studies of the two-leg 
Hubbard model close to half-filling also found that the 
peak of S (q) appears at q=(n, We have also ob- 

served that, as the density decreases from n=l, the peak 
at q=(n, l)n decreases, while another peak at q=(2n, Q)tt 
starts to increase, such that at n ps 0.5 they have the 
same magnitude. In contrast, for large values of J, S (q) 
is one order of magnitude smaller and has a small cusp 
at <f=(2n, 0)7T, as can be seen in Fig.[2td). 

Spin dimerization of the localized spins has been de- 
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Figure 3: (a) Nearest-neighbor spin correlations of the 2-LKL 
for L=32, J=0.8, and n=l/4. Solid and dashed lines repre- 
sent AFM and FM correlations, respectively. The thickness 
of the lines is proportional to the magnitude of the correla- 
tions. Only the ladder central portion is presented. Below the 
correlations, a classical configuration compatible with them 
is shown, (b) Same as (a) but for n=l/2; (c) Dimer order- 
parameter D (L/2 - l)vs. 1/L for J=0.8 and ra=l/4. 



tected in the ID KLM for both J<0l2 and J>0M It 



would be very interesting if the dimerization also survives 
in the 2-LKL, as this would suggest that it may also be 
present in the 2D system. Indeed, we have observed spin 
correlations in the 2-LKL that resemble the dimerization 
of the ID KLM. In Fig.EJa), we show the spin-spin corre- 
lations (S~ • Sf ) for nearest-neighbor sites of the 2-LKL 
at density n=l/4, J=0.8 and L=32. The solid (dashed) 
lines indicate that D(j) is negative (positive) and the line 
thickness is proportional to the amplitude of the corre- 
lations. As can be seen, along the legs the dimer order 

parameter D(j) = (sJ 1 ^ ■ J+1 )^> oscillates with period 

2, while the rungs exhibit FM correlations. We have also 
found that for n=l/2, D(j) also oscillates with period 
2, as shown in Fig.^b). However, in this case, the cor- 
relations along the rungs are antiferromagnetic. This is 
not a finite-size effect artifact or caused by open bound- 
aries: in Fig. die), the order parameter at the center of 
the ladder D(j = L/2 - 1) vs. 1/L, at J=0.8 n=l/4, 
shows a very weak size dependence. For other densities, 
more complex spin structures were observed and there is 
no analogous simple picture (a similar situation occurs in 
the ID KLM away from quarter-fillingifi) . For J^>1 and 

0.5< n<l, (S? - ST) 10- 3 (for n<0.5 some small 

ferromagnetic correlations start to develop), much less 
than the values found for small J. We also have veri- 
fied that, as in the ID case^ these spin correlations can 
be traced back to the RKKY interaction between local- 
ized spins. This effective, conduction-electron mediated 
spin-spin interaction can be obtained from second-order 
perturbation theory and it is given by 



ffRKKY ~ J 



RKKY I 



j|)(Sj • S) + ■ S?) + 2J* KKY (|i - j\)Sl ■ S 



where 



(^rkky(O): ^rkky(I): ^rkky(2), ■••) 



(0, -1, 0.66, 0.07, -0.41, -0.03, 0.28, ...) n = 1/4 
(0,-0.17,1.25,0.73,-1.29,0.10,0.46,...) n = 1/2 



(^rkky(O)j ^rkky(1)> ^rkky(2)j ■■■) — 



(-1.46, 0.23, 0.90, 0.09, -0.41, -0.03, 0.28, ...) n=l/4 
(4.35,2.30,-2.42,0.07,0.70,0.40,-0.86,...) n = 1/2 ' 



r 



and Sj (Sj) is the localized spin in the first (second) leg 
and rung j. For simplicity, only the first few RKKY 
couplings were shown. Let us now focus on density 
n=l/4. All Jrkky(0 nave signs that favor a classi- 
cal configuration along the legs as TTllttll- Note that 
spin dimerization is expected in a spin chain with first- 
and second-neighbor interactions J\ and J2, if J2>0 and 
-4J 2 < Ji<031 The first two RKKY couplings Jr K ky( 1 ) 
and </rkky(2) along the legs satisfy this inequality, and 



further neighbors couplings favor the classical configura- 
tion. Moreover, the couplings Jrkky(0 between legs also 
favor the predominant ferromagnetic alignment across 
the rungs (except for Jj| KKY (l)=0.23, which is neverthe- 
less small), also in agreement with the classical picture 
presented in Fig. Eta) . A similar analysis also holds for 
the case n=l/2. It is interesting to note that, in this case, 
the legs are coupled antiferromagnetically and >/rkky(0 
has larger magnitudes than at n=l/4. This fact suggests 
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Figure 4: The charge structure factor N (q) vs. for the 
2x32 cluster and n=7/8: (a) J=0.8; (b) J=10. We also show 

N<? (<j) and N$ (q) (see text). 



that at n=l/2 the legs are more strongly coupled than 
at 7i=l/4. Indeed, our numerical results of Figs.Efa)-(b) 
confirm this expectation. Thus, the RKKY interaction 
appears to naturally lead to the spin structure shown in 
Figs. Eta)- (b). It is interesting to note that unusual or- 
dered spin structures have been observed in some heavy 
fermion compounds (see, for example, Ref . l2lj) . Our re- 



sults suggest that the effective RKKY interaction may 
be their origin. 

We have also calculated the charge structure 
factor N (q)=^ E '£r 1 ,r 2 e^-^) (5n(f{)Sn(r 2 )) , where 
6n(fi)=n(f[) — (n(f{)}. Previous work on the ID KLM 
has shown that the qualitative behavior of N (q) in the 
extreme weak- and strong-coupling limits could be as- 
cribed to free spin- 5 and spinless fermions, respectively^ 
The same analysis clarifies the behavior of N (q) in the 2- 
LKL. Let us call Nq (q) the charge structure factor of free 
fermions with spin-S" in a two-leg nearest-neighbor tight- 
binding ladder. 2 2. In Fig. Bfa), N (q) is shown for the 2- 

LKL with L=32, J=0.8, and n=7/8, as well 2V Q ' (q). The 
behavior of N (q) is fairly similar to free spin-i fermions. 
On the other hand, in the strong coupling limit N (q) ap- 
proaches the structure factor Nq (q) of spinless fermions 
(see Fig.Etb)). 

In conclusion, we have explored the phase diagram of 
the two-leg Kondo lattice model away from half-filling. 
Our results show that a ferromagnetic phase is present 
only for small densities, unlike the ID Kondo chain but 
consistent with mean field studies^ We have found that 
the charges behave basically as free fermions. On the 
other hand, the spins have non-trivial behavior. The 
peak of the spin structure factor for small values of J 
is located at q=(2n,0)n for n<0.5, and at q=(n, l)ir for 
n>0.5. For large values of J and n>0A S (q) has only 
a small cusp at q=(2n,0)ir. We have also shown that 
dimerization is present in the 2-LKL at densities 71= 1/4 
and n=l/2, and that the RKKY interaction can tenta- 
tively explain this unusual spin arrangement. 
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